arXiv:1509.01527v2 [cond-mat.mes-hall] 8 Feb 2017 


Sound waves induce Volkov-like states, band structure and 
collimation effect in graphene 


M Oliva-Leyva^, Gerardo G Naumis^,^ 

^ Depto. de Fisica-Qmmica, Institute de Fisica, Universidad Nacional Autonoma de Mexico (UNAM). 
Apdo. Postal 20-364, 01000, Mexico D.F. 01000, Mexico 

^ School of Physics Astronomy and Computational Sciences, George Mason University, Fairfax, Virginia 
22030, USA 

E-mail: molivaOfisica.unam.mx, naumisOfisica.unam.mx 


Abstract. We find exact states of graphene quasiparticles under a time-dependent deformation (sound 
wave), whose propagation velocity is smaller than the Fermi velocity. To solve the corresponding effective 
Dirac equation, we adapt the Volkov-like solutions for relativistic fermions in a medium under a plane 
electromagnetic wave. The corresponding electron-deformation quasiparticle spectrum is determined by the 
solutions of a Mathieu equation resulting in band tongues warped in the surface of the Dirac cones. This 
leads to a collimation effect of electron conduction due to strain waves. 



1. Introduction 

Contrary to the parabolic dispersion of charge 
carriers in most materials, the quasiparticles in 
graphene exhibit a linear relation between energy 
and momentum, and thus behave as massless 
relativistic fermions [1]. Consequently, the low- 
energy description for the quasiparticles in graphene 
is given by the massless Dirac equation, with an 
effective “speed of light” of8xl0^m/s. This property 
results in a number of unprecedented features, such as 
Klein tunneling effect, specific integer and fractional 
quantum Hall effect, a “minimum” conductivity of 
^ 4e^//i even when the carrier concentration tends 
to zero, weak antilocalization, high mobilities of up 
to 10®cm^V~^s“^, universal transmittance expressed 
in terms of the fine-structure constant, among others 
[2, 3, 4], 

Graphene, also shows unique mechanical prop¬ 
erties [5]. Suffice it to say that graphene has an 
effective Young modulus of ~ ITPa and simulta¬ 
neously, can reversibly withstand elastic deforma¬ 
tions up to 25% [6]. This unusual interval of elas¬ 
tic response has opened a new opportunity to ex¬ 
plore the strain-induced modifications of the its elec¬ 
trical, chemical and optical properties, and thus, to 
improve its technological functionality. For exam¬ 
ple, a band-gap opening has been achieved by us¬ 
ing an uniaxial strain [7, 8]. On the other hand, via 
stretching of the supporting flexible substrate, pro¬ 
duces impressive increases in the chemical reactivity 
of graphene [9]. Also, the nonlinear response of na¬ 
noelectromechanical graphene resonators has opened 
up new device applications [10, 11]. Very recently, 
the modulation of the transmittance for graphene 
under an arbitrary uniform strain has been quanti¬ 
fied [12, 13]. Moreover, from a view point of basic 
research, strained graphene provides a platform for 
studying exotic properties such as fractal spectrum 
[14, 15], mixed Dirac-Schrodinger behavior [16, 17], 
emergent gravity [18, 19], topological insulator states 
[20, 21], among others [22]. 

Nevertheless, among the most interesting strain- 
induced implications one can cite the experimental 
observation of a spectrum resembling Landau levels 
in strained graphene [23, 24], which was predicted 
earlier by means of gauge fields [25, 26]. Starting 
from a tight-binding elasticity approach, in the 
continuum limit, one can predict the existence of 
these strain-induced gauge fields. Thus, a nonuniform 
deformation of the lattice can be interpreted as a 
pseudomagnetic field, which has been deeply explored 


[27, 28, 29, 30, 31, 32, 33]. However, within 
the same theoretical framework, a time-dependent 
deformation gives rise to a pseudoelectric field but 
its consequences have been less considered [34, 35, 
36, 37]. 

The principal motivation of the present work is 
to determine the consequences of a strain wave (a 
time-dependent deformation) on electron motion in 
graphene. This paper is organized as follows. In 
section 2 we discuss the effective Dirac Hamiltonian 
for graphene under a strain wave whereas in section 
3, we obtain the corresponding solutions. Section 4 
contains discussions and conclusions. 

2. Graphene under a strain wave 

Consider a time-dependent deformation field of the 
graphene lattice u{x,y,t) described by 

It = (0, uo cos(Gy — wt)), (1) 

under the condition of continuum limit where uq <C 
oo 2'!r jG, i.e. the atomic displacement uq is much 
less than the unstrained carbon-carbon distance Oq 
while the wavelength ^-KjG is much greater than 
oq. This deformation wave propagates along the y 
direction with a velocity Vs = oj/G ~ 2 x 10"^m/s, 
which is assumed to be equal to the sound velocity 
in graphene [38]. As illustrated in figure 1 (a), we 
chose the arbitrary coordinate system xy in such a 
way that it is rotated an arbitrary angle 9 respect 
to the crystalline coordinate system xq^o- For the 
latter, the xo-axis points along the zigzag direction 
of graphene sample. Thus, for 0 = 0 the deformation 
wave moves along the armchair direction of graphene 
lattice, whereas for 9 = 7r/2, it moves along the zigzag 
direction. 

As discussed extensively in the literature, the 
electronic consequences of a nonuniform strain can 
be captured by means of an effective gauge field A, 
which in the rotated frame xy is given by [39], 


o 1 

{uxx — Uyy) COS 30 — 2uxy sin 30^, 


o 1 

—2uxy cos 30 — {uxx — Uyy) sin 30^ , 

(2) 


where ~ 3 is the electron Griineisen parameter and 
Uij = {diUj + djUi)/2 is the symmetric strain tensor. 
From the expression for the effective gauge field, it 
is clear that A exhibits a periodicity of 27r/3 in 0, 
which reflects the discrete rotational invariance, i.e. 
the trigonal symmetry, of the underlying honeycomb 
lattice. As is well known, in graphene there are two 
inequivalent Dirac cones, located at special points of 
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Figure 1. (Color online) (a) Schematic 
representation of a sound wave in graphene 
sample. The dark regions represent 
zones of higher density of carbon atoms 
in graphene sample. In the inset, the 
relation between the arbitrary coordinate 
system xy and the crystalline coordinate 
system is shown. (b) Unstrained 

graphene under a pseudoelectromagnetic 
wave: an equivalent description to the 
problem of strain wave. The pseudoelectric 
field lie in graphene plane, whereas the 
pseudomagnetic field is perpendicular to 
the graphene sample. 


high-symmetry of the reciprocal lattice and usually 
denoted by K and K'. The fictitious field (2) has 
opposite signs at different valleys [4] . If for the valley 
K the effective gauge field is A, then for the other 
valley IC' is given by —A. This reflects the fact 
that elastic deformations do not violate the time- 
reversal symmetry [40]. Here, for simplicity we will 
deal with the IC valley, since the other is analogous 
and obtainable by simply changing the sign of A. 

For the case of the deformation field (1), from 
(2), the following effective gauge field results: 

A= sm(Gy — ujt)(cos30, — sm30). (3) 

2ao 

Hence the strain wave (1) gives rise to a pseu¬ 
doelectromagnetic wave propagating along the y-axis 
at the velocity Vg- Here the associated pseudomag¬ 
netic field given by d^Ay — dyA^ oscillates perpen¬ 
dicularly to the graphene sample (see figure 1 (b)). 


On the other hand, the associated pseudoelectric field 
given by —dtA oscillates in the sample plane but, in 
general, it is not perpendicular to the propagation 
direction of the strain wave. Let us illustrate this pe¬ 
culiar behavior. For example, for 0 = ±7r/2 -|- jZ, 
i.e. when the strain wave moves along the zigzag di¬ 
rection, the pseudoelectric field oscillates along the 
propagation direction of the pseudoelectromagnetic 
wave. Thus, for 6 = ±7r/2 -|- 2n'Kl2>, the pseudoelec¬ 
tromagnetic wave (3) behaves more like a sort of lon¬ 
gitudinal mechanical wave. Also, notice that for this 
case the pseudomagnetic field is zero. In contrast, for 
9 = mr/3, i.e. when the strain wave moves along the 
armchair direction, the pseudoelectric fields oscillates 
transversally to the propagation direction of the pseu¬ 
doelectromagnetic wave, as the expected behavior of 
a real electromagnetic wave. 

Then including the strain-induced gauge field via 
minimal coupling, the effective Dirac equation reads 

vpcr ■ {—ihV — A)Af = ihdtAf, (4) 

where vp = 8 x lO^m/s is the Fermi velocity and 
<7’ = (o'xjO'y) are Pauli matrices acting on sublattice 
space. This equation describes low-energy excitations 
of the electronic system in graphene under the strain 
wave. 

3. Volkov-like states 

From (3) one can clearly distinguish the periodicity 
of the effective gauge field A along the time 
direction. Therefore, one could think the standard 
use of the Floquet theory to analyze time-periodic 
Hamiltonians [41]. At the same time, A presents 
space periodicity along the y-direction, so that, Bloch 
(Floquet) theory could also be used. However, since 
A ultimately depends upon the plane phase (j) = 
Gy — Lot of the wave, to solve (4) , we propose a spinor 
wavefunction yp of the form [42] 

^{x,y,t) = eyip[i{k,gX + kyy - Et/K)\^{(j)), (5) 

and then, a Floquet analysis could be translated to 
the effective equation of the spinor, as carried out in 
$(<()) [43]. In relativistic mechanics, this procedure 
is equivalent to a jump into the light-cone frame of 
reference. In principle, E,kx,ky are just parameters 
of the ansatz. However, to recover the solution for 
an electron in undeformed graphene, obtained as 
Ao —>■ 0, it is needed that E = ±hvF{k'/.-\- ky)~^^‘^, as 
is done in similar relativistic problems [43, 44, 45]. 

Our ansatz (5) was firstly used by Volkov to find 
exact solutions of the Dirac equation for relativistic 
fermions under a plane electromagnetic wave in 
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vacuum [42]. It is worth mentioning that only in 
the case when the electromagnetic wave propagates 
in vacuum the solutions of the Dirac equation can be 
found in a simple closed form, these are the Volkov 
states. However, when one consider the interaction 
with a plane wave propagating in a medium with 
an index of refraction ^ 1, the mathematical 
complexity of the Dirac equation is largely increased, 
and in fact, it is a problem with a long history (see 
for example [46, 47] and references therein). 

Comparing our problem, given by (3) and 
(4), with the standard problem of real relativistic 
fermions under a plane electromagnetic wave in a 
medium with an index of refraction rim, one can 
recognize the following. The Fermi velocity vp 
plays the role of the speed of light in vacuum 
c, whereas the strain wave velocity Vs plays the 
role of the electromagnetic wave speed clum in a 
medium. From this analogy, one can think that our 
pseudoelectromagnetic wave propagates in a medium 
with effective index of refraction rim = vp/vg > 1, 
which is a different limiting than the one considered 
in the real electromagnetic case where rim < 1 [48]. 
Only in the hypothetical case that vp = Vg, one can 
obtain the usual Volkov states [49]. 

Substituting (5) into (4) and taking into account 
Vp ^ Vg, we obtain the following differential system 
for the components of spinor $(<()), 

= {kx - iky - sin - E^a, 

dcp 

= {kx + iky - sin()))$A - E^B: (6) 

a<p 

where we define the non-dimensional parameters: 
kx,y = kx,ylG, Aq = /3uo/(2ao) and E = E/{hvpG). 
This system can be reduced to a second-order 
differential equation for each component of spinor 
€>(())), denoted by Fyi and Tb- However, before it 
is appropriate to carry out a new ansatz: 

€>(())) = exp[—iAy())-|-i^osin30coS()>]r(())). (7) 

As a result, from (6) and (7) we obtain that 
the both components F^ and F b satisfy the Mathieu 
equation (see Appendix): 

+ {a-2qcosOT a.b =0, (8) 

where we introduce the variable 
C = (</< +^o)/2, with tan (^0 = 2fca,, (9) 

and the new parameters 

a = Aky, and g = 2Ao(l -I- 4fc^)^/^ cos 36*. (10) 
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Figure 2. (Color online) Stability 

chart as a function of the non-dimensional 
parameters a and q for Mathieu solutions. 
The regions of stability (white domains) 
and instability (red domains) are divided 
by the characteristic curves an{q) (solid 
lines) and bniq) (dashed lines). The chart 
is symmetrical respect to the a-axis. 

Therefore, the general solutions to the compo¬ 
nents r A and T b are linear combinations of the 
Mathieu cosine C(a, q, () and Mathieu sine S{a, q, () 
functions. Nevertheless, taking into account that 
when the pseudoelectromagnetic field is switched out 
(Aq = 0 ) the wavefunction IF must reduce to a free- 
particle wavefunction, we get 

r(C) = Nida, 9, C) + ^s{a, 9 , 0) ( J.a ) , (h) 

where a = tan{ky/kx), V is a normalization 
constant and s = ±1 denotes the conduction and 
valence bands, respectively. To make sure that (11) 
reproduces the case of undeformed graphene, it is 
enough to consider the properties of the Mathieu 
functions for 9 = 0 : (7(0,0, C) = cos(-\/a([) and 
S'(a, 0, C) = sin(^C). 

4. Discussion 

As it is well documented [50] , the stability of Mathieu 
functions depends on the parameters a and 9 . In 
figure 2 , the red regions in the (o, 9 )-plane are those 
for which the solutions are unstable (exponential 
functions) and therefore, are not acceptable from 
a quantum view point. On the other hand, the 
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Figure 3. (Color online) Collimation 
effect of electron conduction by strain- 
waves. (a) Chart of the allowed bands 
(white regions) for the quasi momentum 
k (in units of G), with Aq = 0.15 

and the strain wave propagating along 
the armchair direction. The diagram is 
symmetrical respect to both axes. (b) 
Dirac cone. Red (blue) strips correspond 
to the forbidden (allowed) values of the 
quasi-energy E{kx,ky) (in units of hvpG), 
because of the strain wave. 


white regions are those for which the solutions are 
acceptable wavefunctions. The boundaries between 
these regions are determined by the eigenvalues, 
an{q) and bn(q), corresponding to the 27r-periodic 
Mathieu functions of integer order, cen(q,C) 
cen{q,C), respectively [50]. 

As a consequence from the properties of Mathieu 
solutions mentioned above, a “band structure” 
naturally emerges in our problem. If one take into 
account (10), the stability chart in the (a, ( 7 )-plane is 
translated into a chart of the allowed bands in the 
{kj;,ky)-p\a,ne, as shown in figure 3 (a). When the 
strain wave propagates along the zigzag, cos 30 = 0, 
then q = 0 and therefore the quasi-wave vector 
fc = ikx,ky) can take any value. However, from the 
dependence q ~ cos 30, one can conclude that for a 
strain wave propagating along the armchair direction 
(cos 30 = ±1), the band gaps in the {kx, fcy)-plane are 
expanded. Just for the sake of illustration, hereafter 
we consider 0 = mr/S. 

In figure 3 (a), we display the allowed values 
of the quasi-wave vector k = {kx,ky) for a strain 
wave of amplitude uq = O.Ioq and moves along 
the armchair direction. The most important result 
to notice in figure 3 (a) is that the quasi-particles 
propagate preferably in the y direction, i.e., in the 
propagation direction of the strain wave leading to 
a collimation effect of the electrons. For example, 
notice that if ky = 0, the allowed values of kx are 


practically limited to the interval (—gc/4Ao, 9c/4Ao), 
where bi{qc « 0.91) = 0. In other words, only low 
energy quasi-particles propagate perpendicularly to 
the propagation direction of the strain wave. On the 
contrary, for kx = 0, \ky\ can take all values except 
basically those of the form to/ 2, where to is a positive 
integer. This last result is an analogous condition as 
Bragg’s diffraction. 

To end, let us point out that the forbidden values 
of the quasi-wave vector k divide the Dirac cone, 
given hy E = iQtx + kyY^'^, into strips of allowed and 
forbidden values of the quasi-energy E, as illustrated 
in figure 3 (b). Similar results for the band structure 
have been discussed in earlies works for the case of 
a spin-less particle (which obeys the Klein-Gordon 
equation) in a medium (of Um > 1) irradiated with 
an electromagnetic plane wave [43, 44]. On the 
other hand, our results differ from those reported for 
graphene under an electromagnetic plane wave moves 
with the velocity of light in vacuum c » 300?;^ [48]. 
The physical reason is simple. In our problem, the 
pseudoelectromagnetic wave moves with the velocity 
of sound Vs vp/AO. 

Vaezi et al [36] reported the possibility of 
observing a nonvanishing charge current in graphene 
by applying a time-dependent strain. To archive this 
edge charge current, it was assumed a mass term 
to provide the presence of a gap in the spectrum, 
which is an essential ingredient to get a quantized 
response. From figure 3 (b) one can be distinguished 
that graphene remains gapless. Therefore, it does not 
seem possible to archive such edge charge currents in 
graphene under the considered approximations. To 
achieve a gap and thus topological modes as happens 
in the electromagnetic case [48], this will require a 
much higher amplitude of the strain than the one 
considered here since intervalley mixing is needed. 

This is a consequence that in fact, there is a cru¬ 
cial difference between magnetic and pseudomagnetic 
fields. As commented above, the coupling constant 
for pseudomagnetic fields has opposite signs for elec¬ 
trons at different valleys, whereas for a magnetic field 
is the same. Since the coupling constant is valley 
anti-symmetric, the currents at the K and K val¬ 
leys flow in the opposite directions and they cancel 
out [37]. As a result, in general a pseudoelectric field 
does not cause a net electric current and thus adia¬ 
batic pumping is not possible as happens in the real 
electromagnetic fields [48]. Notice that in our case, 
the pseudomagnetic field can be considered as adia¬ 
batic since w <C Gvp- 
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5. Conclusion 

In conclusion, we studied the effects of a strain 
wave on electron motion in graphene. The 
coupling between the quasi-particles and the strain 
wave are captured by means of an effective 
pseudoelectromagnetic wave. As solutions to the 
resulting effective Dirac equation, we found Volkov- 
type states, which propagates preferably in the 
propagation direction of deformation. Also, we 
reported a band structure of allowed and not allowed 
values for the quasi-momentum and for the quasi¬ 
energy. The form of the emergent band structure 
depends on the propagation direction of the strain 
wave respect to the crystalline directions of graphene 
lattice. This fact produces a collimation effect 
of charge carriers by strain waves, which should 
be an alternative mechanism to archive electron 
beam collimation, beyond magnetic focusing [51], 
an external periodic potential [52] or nanostructured 
heterodimensional graphene junctions [53]. 
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Appendix 

In this section, we present the details of the 
calculations to derive (8) from (6). Note that, the 
differential system (6) can be rewritten as 


+ ~ ^0 


I -I- 4 fc 2 cos 30 COs{(j) + (/>o)]rA,B = 0. 


d(j) 


= {kx — Aq cos 30 sin 0 — iky — iAo sin30sin((i)<i)B 




d^A 

dcj) 


= {kx — Aq cos 30 sin 4> + iky -I- iAg sin 30 sin (j))^A 


—E^b- 

To simplify this system, one can propose that 
f4’ 

= exp[ / (—— *Aosin30sin0*)d(/)*]r(^), 

= exp[—iky(j) + iAo sin 30 cos (j)]r{<j)), 
and then one get the following differential system 

dr B 


d(j) 
dVA 
d(j) 


= {kx — Aq cos30sm(p)TB — ETa, 

= {kx — Aq cos 30 sin 4>WA — EVb, 


for the components of the spinor r(^). Now, taking 
second derivate one can reduce this last system to 
a Hill differential equation for both components P^i 
and P^: 

— [1’^ -I- [A^ — Aq cos 30 cos (j) — {kx — Aq cos 30sin0)^]r^ s = 0. 
d(p^ 

However, since in the our problem Aq = 

/3uo/{2ao) ^ I, the second-order terms in Aq can be 
neglected, and taking into account that ky = E^ — kx, 
then one find that 
d^TA^B 
d(jA 

where tan^g = “^kx- Finally, if one carry out the 
variable change ^ = ((/> -I- <(>o)/2 and define the 
parameters a = Aky and q = 2Ao(l -|- cos 30, 

one obtain (8), which is the Mathieu equation. 
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